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PREFACE 


This  paper  is  the  ninth  in  a series  of  reports  growing  oat  of  studies 
of  radar  cross -sections  at  the  University  of  Michigan’s  Willow  Run 
Research  Center.  The  primary  aims  of  this  program  are: 

1.  To  show  that  radar  cross-sections  can  be  determined  analyti- 
cally. 

2.  To  elaborate  means  for  computing  cross-sections  of  objects  of 
military  interest. 

3.  To  demonstrate  that  these  theoretical  cross-sections  are  in 
agreement  with  experimentally  determined  values. 

Intermediate  objectives  are: 

1.  To  compute  the  exact  theoretical  cross  -sections  of  various 
simple  bodies  by  solution  of  the  appropriate  boundary-vaiue 
problems  arising  from  the  electromagnetic  vector  wave 
equation. 

2.  To  examine  the  various  approximations  possible  in  this  prob- 
lem, and  determine  the  limits  of  their  validity  and  utility. 

3.  To  find  means  of  combining  the  simple -body  solutions  in  order 
to  determine  the  cross-sections  of  composite  bodies. 

4.  To  tabulate  various  formulas  and  functions  necessary  to  enable 
such  computations  to  be  done  quickly  for  arbitrary  objects. 

5.  To  collect,  summarize,  and  evaluate  existing  experimental 
data. 

Titles  of  the  papers  already  published  cr  presently  in  process  of  publi- 
cation are  listed  on  the  back  of  the  title  page. 
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NOMENCLATURE 


A 


= semi-major  axis  of  the  oblate  spheroid. 
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= a condensation  symbol  --  Equation  (44). 

= semi -minor  axis  of  the  oblate  spheroid. 

= condensation  symbols  — Equation  (77). 

= Beta  function. 

- the  oblate  spheroidal  coefficient. 

--  incident  electric  field  vector. 

= x,y,z  components  of  the  incident  electric 
field  vector. 

= the  incident  electric  vector  E*  with  the  phase 
iscior  removed. 


E,  O 
F 


= magnitude  of  the  incident  field  vector  E\ 

= denotes  whether .i?  and  L are  even  or  odd  in 
the  In(^,L). 

= 1/2  the  focal  length  of  the  coordinate  ellipsoid 
and  hyperboloid. 

= incident  magnetic  field  vector. 


H*  .H1 
-x  — y 


= x,y,z  components  of  the  incident  magnetic 
field  vector  2*. 


H 


C 


- magnitude  of  the  incident  magnetic  field 
I 

vector  H . 
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H, 


lnU'.U 
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m 
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R(n),S(0,#(*) 


VJm«> 


= the  incident  magnetic  field  vector  H*.  with 
the  phase  factor  removed. 

= integrals  occurring  in  Equations  (68)  and 

(69). 

= a solution  to  the  vector  wave  equation  con- 
structed from  solution  to  the  scalar  Helm- 
holtz equation. 

= the  norms  of  the  angular  functions. 

= a point  of  observation  in  space  given  in 
either  oblate  spheroidal  coordinates  (>?,  £,<*>) 
or  spherical  coordinates  (r.0,<*>). 

- total  power  intercepted  by  an  isotropic  scat- 
terer  of  cross-section  cr. 

= Foynting  vector  of  incident  electromagnetic 
field. 

= Poynting  vector  of  scattered  electromagnetic 
field. 

the  location  of  a radiating  source  in  space 
given  in  either  oblate  spheroidal  coordinates 
(*?' .£'»£')  or  spherical  coordinates  (r ' ,©*,<#>*). 

= solutions  of  the  separated  scalar  Helmholtz 
equation. 

= angular  oblate  spheroidal  wave  function. 

= radial  oblate  spheroidal  wave  function. 
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NOMENCLATURE  (Continued) 


= an  index  used  to  indicate  any  of  the  x,  y,  or 
z coordinates . 

- n vector  used  to  indicate  any  of  the  unit 

vector  s i , i , or  i . 

-x  ~y  ~z 

= the  coefficients  in  the  power  series  expansion 
of  r, 


= the  imaginary  unit  V - 1 . 


A -i,  >1 
x y z 


unit  vectors  in  rectangular  Cartesian  co- 
ordinates . 


-ri  ’h  ‘-4> 


unit  vectors  in  oblate  spneroidal  coordinate 
system. 


= r :r.  tnth  rrr» r* 


i_  /\ 
*-11/  A , 


a separation  constant  which  is  found  to  be  an 
integer . 


^u(c  {) 


= the  spherical  Neumann  function. 


a ,i 
Am 


a/'P! 


= the  coefficients  in  the  expansion  of  the 
radial  functions  of  the  second  kind. 

= the  radial  distance  from  the  origin  to  a point 
in  space. 

- the  differential  operator  i i + i . 

_x  <)x  ~y  8 y ~z  dz 

- arbitrary  constant  determined  by  boundary 
conditions . 
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NOMENCLATURE  (Continued) 


Vm 

ylm 


n 

\ 

* 

p 

c r 


4> 


- a separation  constant. 


= a separation  constant. 

= the  Kronecker  delta. 

= kF. 

= coordinate  hyperboloid  of  the  oblate  sphe- 
roidal coordinate  system. 

- wavelength. 

= coordinate  ellipsoid  of  the  oblate  spheroidal 
coordinate  system. 

= distance  between  two  points  P and  Q in  space. 
= scattering  cross-section. 

- coordinate  angle  of  the  oblate  spheroidal 
coordinate  system  measured  counter- 
clockwise around  the  z-axis  from  the 

xz -plane . 

= a solution  of  the  scalar  Helmholtz  equation. 
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INTRODUCTION  AND  SUMMARY 


Ac  the  present  time , exact  solutions  are  known  for  the  electro- 
magnetic scattering  cross-sections  of  only  four  three-dimensional 
configurations:  the  sphere,  the  paraboloid,  the  semi -infinite  cone, 
and  the  prolate  spheroid.  In  this  paper,  the  exact  back-scattering 
cross-section  of  a fifth  three-dimensional  configuration,  the  oblate 
spheroid,  is  obtained.  Although  this  body  may  be  of  little  practical 
importance,  the  scattering  cross -section  has  been  obtained  as  part 
of  the  first  intermediate  objective  stated  in  the  Preface,  to  widen  the 
theoretical  knowledge  of  this  field.  The  solution  is  obtained  for  one 
direction  of  incidence  and  presented  in  the  form  of  an  infinite  series. 
Much  work  still  remains  to  be  done  in  developing  methods  for  evalu- 
ating this  series  to  obtain  numerical  values  of  cross-section. 

The  electromagnetic  scattering  cross-section  of  the  oblate  sphe- 
roid has  been  obtained  as  follows:  the  oblate  spheroidal  coordinate 
system,  in  which  the  scalar  wave  equation  is  separable  and  which 
includes  the  scattering  surface  as  a coordinate  surface,  is  chosen. 

A plane,  linearly  polarized  wave  progressing  along  the  z-axis  in  the 
negative  z -direction  is  incident  on  the  surface  of  a perfectly  conduct- 
ing oblate  spheroid.  Using  the  method  of  Hansen  in  which  a solution 
of  the  vector  wave  equation  is  constructed  from  the  solution  of  a scalar 
Helmholtz  equation,  an  expression  is  obtained  for  the  scattered  electric 
field  in  series  form  having  a number  cf  arbitrary  constants.  This  ex- 
pression is  valid  everywhere  in  space.  By  means  of  the  boundary  con- 
ditions at  the  surface  of  the  scatter er,  the  properties  of  the  scattered 
radiation  are  related  to  those  of  the  incident  radiation  so  as  to  obtain 
defining  relations  for  the  arbitrary  constants.  Unfortunately,  these 
defining  relations  are  quite  complex  and  no  simple  technique  has  vet 
been  developed  for  determining  the  exact  values  of  the  required  con- 
stants. 

The  scattered  electric  field  is  used  to  determine  the  hack-scat- 
tering cross-section  of  the  oblate  spheroid.  The  back-scattering 
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cross-section,  c.  is  defined  by  Ridenour  (Ref.  i)  as:  V (dimensions 
of  an  area)  is  to  be  4tt  times  the  ratio  of  the  power  per  unit  solid 
angle  scattered  back  toward  the  transmitter  to  the  power  density 
(power  per  unit  area)  in  the  wave  incident  on  the  target.  In  other 
words,  ii  at  the  target  the  power  incident  on  an  area  cr  placed  nor- 
mal to  the  beam  were  to  be  scattered  uniformly  in  all  directions,  the 
intensity  of  the  signal  received  back  at  the  radar  set  would  be  just  what 
it  io  in  the  case  of  the  actual  target.1'  The  expression  obtained  for  the 
scattering  cross-section  depends  upon  the  set  of  oblate  spheroidal  co- 
efficients which  have  not  been  tabulated  extensively  and  upon  the  above- 
mentioned  arbitrary  constants.  Therefore,  the  usefulness  of  the  ex- 
pression for  the  cross-section,  as  far  as  numerical  calculations  are 
concerned,  is  limited  at  the  present  time. 

In  conclusion,  the  author  wishes  to  acknowledge  his  indebtedness 
to  the  work  of  F.  V,  Schultz  (Ref.  2)  and  of  A.  L,eitner  and  R.  D.  Spence 
(Ref.  3). 
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II 

THE  OBLATE  SPHEROIDAL  FUNC  TIONS 

To  obtain  expressions  for  the  incident  and  scattered  electro- 
magnetic field  and  hence  the  back-scattering  cross-section,  the 
properties  of  the  solutions  of  the  vector  and  scalar  wave  equations 
for  the  spheroid  are  required.  These  properties  are  presented 
below. 

2.  1 THE  COORDINATE  SYSTEM 

The  solutions  of  the  vector  wave  equation  and  of  the  corresponding 
scalar  wave  equation  can  be  expressed  in  oblate  spheroidal  coordinates 
rj , and  </.„  This  system  consists  of  three  families  of  orthogonal  sur- 
faces which  are  described  by  vt  t and  equal  to  constants.  The  first 
two  of  these  families  consist  of  confocal  surfaces  of  revolution  about 
the  z-axis,  namely,  hyperboloids  of  one  sheet  and  oblate  spheroids  as 
indicated  in  Figure  1 . The  equation  <t>  equals  a constant  represents 
the  family  of  half  planes  through  the  z-axis,  where  <f>  is  the  angle  be  - 
t'.vesn  these  planes  and  the  xz*  plane. 


Oblate  Spheroidal  Coordinates 
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The  equations  of  transformation  from  the  oblate  spheroidal  co- 
ordinate systems  to  the  Cartesian  svstem  are: 


x = f| 

[«. 

2 2 ] 

i 1/2 

- *5  M * i-  ? ) j 

j cos  P. 

< 

ic> 

i 

Z 2 1 

1/2 

y — F j 

! 

f 

l«N 

l 

sin  d. 

z = F-je, 

where  F is  1/2  the  focal  length  (PF;),  -1<  *j<1,  arid  0<^<2nr. 

In  the  Cartesian  system  the  equation  of  the  oblate  spheroid  f = f is 

0 

2 2 2 
x + v z 

.2  + ti2  = *’ 

where  A and  B are  the  semi -major  and  semi-minor  axes,  respectively, 
ana  are  given  by 


a-fK!  + «v  b*fv 

The  distance  from  the  origin  to  any  point  in  soace  is  specified  by 

in 


- 

i +•  * - 


1/  a c c | s 

r = I x y + z - F 

When  ? is  large,  r is  approximately  F$. 

The  relations  between  unit  vectors  in  the  Cartesian  system  and  the 
oblate  spheroidal  system  are 

f 212  2 2 -1  2 ] 

1^=1^  {_-*(!  - I ) (?  + «?  ) ] COSO 

r 2 1/2  2 2 -1/2 1 

(**  + i~)  j cos  ^ - j.  sin#  (2a) 

- f ,,  2 1/2,  2 2 -l/zl  . 

ly  = [-*$(1  + e ) (v  r * } 1 j sin  © 

. f „ 2 1/2  , 2 2-1/2 1 . , , , 

+ - *s  } in  + * ) ] snao  + cos  (2b) 
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. ' ,,  2.1/2.  2 2-1/2]  . f ,,  2,1/2,  2 2 -1/2  “I  , , 

-z  = -n  " Tl  ’ (v  + * ) J + h L^1  + ^ f 4 ) J (2c) 


and 


. r 2,1/2  , 2 2-1/2  1 

2r;=ixj'??fi  + 0 + J ) cos  0 J 

. f ,,  2,1/2  - 2 2, -1/2  1 

f I j - j;  (I  + $ ) ' (’i  + f ) sir.  0 ! 


. r * /,  2,1/2 . ?.  2, -i/i 

iz  M(!  - •»  ) U?  + « ) 


r , , 2 1/2  2 2 -1/2  1 
l = l [ ? (1  - n ) ' (r?  + { ) cos  ^ J 


. I"  ,,  2,1/2  , 2 2. -1/2  . 1 

+ i i(l  - r,  ) ( r)  + $ ) sin  <i>  j 

y l j 

F ,,  2 1/2  , 2 2-1/2 

+ iTT|_r'(1+^)  ('J  + * ) 


(3a) 


(3b) 


l , = -i  sin  0 + i cos  0 

-0  -x  ~y 


(3 


* r*  \ 


2.  2 THE  SOLUTIONS  OF  THE  SCALAR  AND  VECTOR  WAVE  EQUATIONS 


The  scalar  Helmholtz  equation. 

\7  '■!■'  + =0, 

in  oblate  spheroidal  coordinates  is 


& 

8tj 


, 2.5'/]  / 1 

+ * ’aTJTT T 


(4) 


1 \ 3 * .2,  2 2, 

7 j — 7 + f W t f )0  = 0 

1 + t;  ) d* 


P.tr  e ^0 

where  * = kF  = — -F,  k = „ , and  X denotes  wavelength 


(5) 


Equation  (5)  can  be  separated  by  assuming  a solution  of  the  form 

tin.  t.d>)  = H('l)  S({)  s (0).  (6) 
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The  substitution  of  Equation  (6)  into  (5)  yields 

r 2 


du 

(1 

d 

— 

d« 

(1 

d~0  (0) 

' ' d 77 

2 dS(t) 
5 >-^7 


» • 


2n 


L'-n 

r 2 


- Q+f  (l-rj)jK(r))=0 


’J 


(?) 


m 2 2 

| T | J - a + t vl  + l } 

J 1 


+ m 0(0)  = 0 


S (5)  = o (8) 


(9) 


where  m and  a are  introduced  as  constants  of  separation.  Since  it  will 

be  convenient  later  to  express  a in  terms  of  m and  another  constant , XL , 

a is  written  as  a . . Particular  solutions  of  Equations  (7)  and  (8)  will 

-f  rn 

be  denoted  by  R{t})  = XT  * and  SU)  = Vym($).  The  functions 

and  V are  known  as  the  angular  and  radial  functions  respectively. 

A solution  of  Equation  (9)  is  0(0)  = cos  m0.  To  insure  single  valuedness, 
m is  restricted  to  integral  values. 

Hansen  {Ref.  4)  has  shown  how  a solution  of  the  vector  wave  equation 


VV  • V - V X (\7x  V)  + k V = 0 


(10) 


can  be  constructed  from  a solution  to  the  scalar  Helmholtz  equation  by 
using  the  solenoidal  relation 


M(a)  - 

iVl  = 


X a 0 1!>  X a 


(ID 


where  a^  is  an  arbitrary  constant  vector  and  4>  is  a solution  to  the  scalar 
wave  equation.  Furthermore,  the  vector  V can  be  expressed  in  terms 
/a\ 

of  M'  . In  particular,  a is  restricted  to  the  Cartesian  unit  vectors 
i , snd  i^  and,  correspondingly,  Mv  ' is  restricted  to  M , M , and  M . 
The  function  ^ is  a Lan'i  product  of  the  functions  R (q ) , S (.*),  and  #(0), 


I 

1 

j 
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2 . 3 EXPRESSIONS  FOR  THE  ANGULAR  FUNCTIONS 

Leii.ner  and  Spence  (Ref,  3)  have  presented  a thorough  discussion  of 
the  properties  of  the  angular  functions.  Since  these  properties  are  used 
quite  extensively  here,  they  are  quoted  verbatim. 

“Equations  (7)  and  (8)  may  be  obtained  from  the  general  equation 


a r„  - 


, d F / ■ ~] 


k.  m ! 


I - h 2 -Vn,  + ' (1  ■'  1 Flm-°  <12> 

J L*  “ z J 


by  replacing  the  complex  variable  z with  fj  or  i respectively.  Equa- 
tion  ^12)  differs  from  trie  sc^u&tion  of  t Vie  sissocicited  H.c^endre  functions 
only  in  that  it  possesses  an  irregular  singular  point  at  z = «.  The  ex- 
ponents associated  with  the  regular  singular  points  at  z = +1  are  fm/2. 

A solution  corresponding  to  the  positive  exponent  will  be  called  a func- 
tion of  the  first  kind:  the  other  independent  solution  will  be  called  a 
function  of  the  second  kind.  Since  the  exponent  difference  is  a.n  integer 
the  solutions  of  the  second  kind  will  possess  logarithmic  singularities  at 
z = ji,  Since  this  point  corresponds  to  ^ = ±i  or • S = + i it  lies  in  the 
physical  range  of  t*  but  outside  the  physical  range  of  ? . Therefore,  the 
angular  functions  are  functions  of  the  first  kind  while  the  radial  functions 


f i-  - ./■  i — i u : i r: a. j j u: 3 

die  iunuuuua  ui  uum.  me  ui  at  emu  octunu  muu* 


’’The  analogy  between  (7)  and  the  equation  for  the  associated  Legendre 
function  suggests  series  of  the  form 


„ , , 2. m/2  > _/m 

- 'I  > A C2k  (1 

k = 0 


2.k  , 

; (.£  - m)  even. 


U/m(^  = M1  -!\2>m/2  2^  Cfkm  (l-^r  KJ-  ~ rn)  odd. 

k = 0 


Also  by  analogy 


. iU  + ym 


I 
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where  v , is  a function  of  e . As  in  the  case  of  the  associated  Legendre 

i.  m 

equation,  Ji  is  restricted  to  integral  values  equal  to  or  greater  than  m. 
This  notation  follows  that  of  Page  (Ref.  5).  Stratton  et  al.  (Ref.  6)  use 
the  symbol,  x , to  denote  (.£  - m)  as  used  here.  The  recursion  formulas 
for  the  coefficients  in  (13)  and  (14)  are 


2k  (2k  + 2m)  C 


k 


m 


2k 


(2k  + m - 2)  (2k  + J+  m - 1)  - r 1 

jt  m j c 


Jm 

2k-2 


2 ^/m 
2k-4 


- t C*  . = 0,  (Jl-  m)  even,  (16) 


’k (2k  + 2m)  - ^(?k  - / + m - 1)  (2k  4 / + m)  - rJm~\  C 


J'm 

2k -2 


2 m 
'2k-4 


cr."\  = o,  (/  - m)  odd  (17) 


with  C = 1 in  either  case*1  The  magnitude  of  the  C J begins  to  de- 

u K 

crease  when  2k>(  x - m)  + 2 if  (X  - m)  is  even  and  for  2k  >(X  - m)  + 1 if 
(X  - m)  is  odd.  For  large  values  of  k 

~1TL 


c 


2 k 
2k-2 


1 


+ /cJm  \ 

2(  C2k-2  \ 
4k  a I 


(18) 


k-2 
Iva 
'2k-4 


!These  recurrence  relations  were  omitted  by  Lcitner  and  Spence. 

Their  continued  fraction  expansions  (Eqs„  (19)  and  (20))  have  been  corrected. 


4(1  + mj  C 


- jVn  -^)(m  + /+  1)  - >•/»  1 cfm  = 0,  ( li-  m) even. 


2 l 
Urn 


4(1  + m)  c‘m  - I <m  - f + !>!m  + J+1)-  rjm  j Cglm  = 0, 

(/- 


m)  odd. 
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regardless  of  - m).  As  k approaches  infinity  this  ratio  must  either 

2 2 

approach  unity  or  approach  zero  as  -«  /4k  . In  the  latter  ease  the 

x-  m 

C_,  alternate  in  sign  at  large  k and  the  series  (13)  and  (14)  converge 

2 K 

absolutely  for  all  finite  r,.  From  (lb)  and  (1?)  one  may  easily  obtain  the 
continued  fractions 


y m = - (i'+m)ti-m) 


/m  " u f , 7 

■ Z{2£  -1}  + y.  + (^+m-2)(/-m-2)  f , 

^ m 4(2^ -3)  + y^  + (/+ m-4)(^-m-4) 

6(2^ -5 )+  y»  + 


+ (j'+m+2)[J-Tn+2)  t ^ 

2(2r  +3)  - y - + 4)(i'-mt  4)  e 

4(2^ +5)  - v a - (Am+6)U'-m-t-6)  « 

^ m " 6(2i+7)  - y/;  + 

-*■  m • 


C /-m)  even. 


ti'+m-l)Ci;-m-l) « 7 

"2(^-D  + >>_  + C-/-rrn-3)(/-m-3)  t" 

4(2y-3)  + y n + (/+m-5)(./-rn-5) « 

im  6<^-5»+Vm 


+ C/+n-.+l)(J?-m+l) 


2(2^+3)  - + (Am+3)(/-m+3)  * 


+ (^+m+5)(j;7-m+5) 
6(2^7)  JTj— 


(Jl-m)  odd.  (20) 
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In  both  (19)  and  (20)  the  first  continued  fraction  contains  only  a finite 
number  of  terms  while  the  second  contains  an  infinite  number  of  terms. 
Equations  (19)  and  (20)  constitute  transcendental  equations  whose  roots 
are  the  eigenvalues  of  >■  .j  . The  eigenvalues  can  be  computed  by 

cutting  off  the  infinitely  extending  continued  fraction  and  solving  the 
resulting  algebraic  equation  by  successive  approximations.  By  re- 
taining more  and  more  terms  in  the  continued  fraction  one  can  obtain 
numerical  values  for  yo  -c  any  desired  degree  o-’  accuracy.1 

2 

'‘A  scries  for  the  v > in  oowers  of  e “ may  be  obtained  (Ref.  5) 

, . ...  ,. 

by  substituting 


Y Jm  2k 
A f2k  • 

k - 0 


into  (19)  or  (20)  and  collecting  like  powers  of  e . One  finds 


,,'m  1)  +■  m - 1 

2 = “ (22  - l)(2Y  + ir 

Pm  - 1 ) H (J)  HU't  l)Il(,^+2) 


2(2.kf  - 1} 


2(2 ,J+  3) 


4m'  - 1 j"  H lj-  1 )(H  (J)  H(/-r  1)  ti(J+  2)  ~| 


r-  2 2 f 

Jm  _ H(/-  1 )H(  4 4(4rn  - 1) 4 

8 2( 2js  - l ) j2^.;-5)2(2^-  l)4  (i^+3)2  2(22.-l) 


H(i-3}HU  - 2) 

+'  2-4(a/-  l)(2^-3) 


'-"A  method  for  rapidly  finding  the  roots  of  such  continued  fractions 
was  independently  developed  by  Bouwkamp  (7)  and  Blsnck  (8).” 


; 


i 

1 


r 

i 


1 

i 


W I LEOVC'  RUN  RESEARCH  CENTER  - UNIVERSITY  OF  MICH 


HCM)Hjj+Z) 
2(2//+  3) 


m 


4(4m 


1) 


[lj - 1 )2  (2y+  3 )4  ( 2/+  7) 2 


H(i+3)H(i+4) 

+ 2- 4(2/+3)(2i+^) 


+ 


2(2/- 3} 


wiigrs 


H(x)  = 


2 2 
x - rr. 


4x  - 1 


(25) 


(26) 


The  series  given  in  (22)  docs  not  appear  to  converge  sufficiently  rapidly 
to  be  useful  if  t is  greater  than  5.  For  smaller  values  of  « it  may  be 
used  to  obtain  a good  approximation  for  y ••  which  may  then  be  im- 
proved  by  use  of  (19)  and  (20). 


" the  angular  functions  corresponding  to  a given  e are  orthogonal  in 
the  interval  (-1,1).  that  is 


,1 

! u; 

^ - ; 


6 , N 

j L>  jam 


(27) 


where 


V im  im  ?m+1+J+1  - H-iV 

N -i  = ? / C CT  - ‘ + / i -rr ) 

im  A*  2i  U2j  l-3-5***2(m  + if))  + l’ 


even . 


I --  O j - o 


OO  no 


(?Z\ 

* t 


i'm  im  2“'‘,J,‘  (m  + i + j)i  , //  % 

^ ■ C->‘  C TTTfe  1' . Tjs 7—  TTT7T7~  , U -m*  odd . 


i = 0 j = 0 


2i 


2j  1 • 3 ' 5 * • *|2(  m + i + j ) + 3j 


(2°) 


2.4  EXPRESSIONS  FOR  THE  RADIAL  FUNCTIONS 

Both  independent  solutions  of  the  radial  equation  (H)  must  be  in- 
corporated ir.  the  general  solution.  For  the  functions  of  the  first  and 
second  kind  Leitner  and  Spence  give 


I GAN 


i 


i 


I 


I 
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<»..  ...  _ (i  +r)m/2  r j. 

v tt~.  / c2 


m 


k - 0 


m 

k 


2 k 

(1+4  ) , even,  (30) 


m 

' *V,7  (*) 


. , 2.ni/2 

til  + { ) ' 


— - tt  /i  /n> 
d rj  ” ' 


v 

vr  = n 


— X. 

^ 2k 


{l  + £2)k,  (/-m)  odd,  (31) 


and 


(2) 


'v  (5)  = 


.i  m 


,,  ,2ni/2  V"  „m+k  , , ^ An  m+k. 

Vm(1  + £ > 2 <m+k)l  CC 

k = 0 ( t£, ) 


U%) 

r 

(-2-m)  even,  (32) 


(2) 


n • •(<?) 


A/ 


, ..  , 2 m/2  * m+k  _ /m  m+k+1 

;Jf)  - »Vm(l  + « 1 * ’ C*  “T75« 

K — U 


Uf) 

(x-m ) odd. 


m+k 


(33) 


where  n represents  the  soherical  Neumann  function,  and  the  con- 

stants  q are  given  by 

,m+k,  . ..  „Jlm.  -1 


q/m 


f ^ 2m+k(m+k)l  C*m  'j 

= l A \ 

■ ^-o  1.3*5...  [2(m+k)+ll  [ 


^ “ '“iiSV"  ^2k  { 


r 


q/m  *j  1*  3*5**»[2(m+k)+3]  ^ 


(A-m)  even. 


(^’-m)  odd. 


* O A \ 


(35) 


Ir.  the  problem  being  considered,  the  primary  interest  is  in  a combination 
of  those  solutions  which  will  represent  a diverging  wave  at  large  distances 

(3) 

from  the  spheroid.  Leitner  and  Spence  show  that  V a (f) . given  by 

(2)  HI  . t z> . . Jirr> 

Vj  (5)  - ' ,V^>^(£)  + i v represents  such  a diverging  wave. 


12 


i 

s 

I 

I 
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They  also  give  an  asymptotic  expression  for  this  function  which  will  be 
very  useful: 


m , ,,m+l  e* 


lim  ' 'V  • (S)  - q y (Z/O  (-i) 

*5>oo  - m ' 


V />  tr\  = r,  . /2/c)m{-i)rn+2  ~ 

7m‘>'  7V"  ' 


■ , even. 


/ Z'  — \ 3 -i 

\^x-  - ha  ; uuu  , 


(36a) 


i 

oou ; 


cl  (3),,.  , , , ,rn  , , ..m  e ? , a , ,, 

lim  — — V/  (J)  = a i2 /« } ml  (-x)  , (£-m)  odd  or  e^en. 

£*->co  7 m 7 £ m f5 


1 

) 
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III 


THE  INCIDENT  ELECTROMAGNETIC  FIELD 


I 


Geometry  Used  in  Describing  the 
Incident  Electromagnetic  Field. 


Let  Q be  a source  point 
and  P a point  of  observa- 
tion given  by  oblate  sphe- 
roidal and  spherical  co- 
ordinates, P or 

Q {f\.  ^ > or 

r',9'.-0').  If  p is  the  dis- 
tance between  P and  Q 
and  y the  angle  between 
r and  r' , then 


P 


1 

r~ 


+ (r') 


2r  r'  cosy 


11/2 

J ' 


(38) 


and 


cos  / = sin  9 sin  O'  cos(0  - o')  + cos  9 cos  6*.  (39) 

If  the  source  point  Q is  moved  to  infinity  along  a fixed  direction  i/, 
0',  then  for  a plane  wave  progressing  aTong  a specified  direction  (Ref.  3 


v 1 2(-i r 


-ikrcosr  \ T i Z(-lp  ..  , v TT 

e = ./  ; t O U !}„  (h)  b /?„  (h!)  V/>n{£,) 


JS,  .J  A 

JL.  in  rn 

> _2  (-1)  i < 


rr0  -*u  ' **  a 

y zi 

^ %%  mi 

'i 

X VJm  <\>  Uim  «V’  '‘’VJ m <F>>  cos  - *•)] 


m : 


(40) 


where  the  N n are  given  by  Equations  (28)  and  (29). 
a.  m 


I 


14 


! 

! 


■i 

i 


I 


GAN 


! 


<4 


j 
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The  problem  is  simplified  if  ine  source  noinr  is  moved  to  infinity 
rlcng  the  positive  z-axis  in  which  caserj  ! = + 1,  z = r cos  r and  Et^aa- 
tion  (40)  becomes 


J 


V 2(-l)“  ,r  , . (1)„ 

N/i(Vo  A ' /ro  k' 


(41) 


where 


JN 


Jn  /n  2i+j+1{i-n>' 

tk  <k  c*  <V  rate  ■ u-°> even-  <42> 


and 


Kio 


Let 


i = 0 j - 0 


A ci°  c{°  ? j AJ  UV  * U-0)  odd.  (43) 

fTo  fTo  21  2J  l^—[2u+])+3j 

J 


2(-l) 


“j/ 0 N /i  n i 

10  *J^0 


and 


(1)Wrs)  - u/o 't1  <1>vio (s) 

i then  Equation  (41)  may  be  written  in  the  form 


(44) 


(45) 


-1KZ 


(1) 


AJo ' Vo 


(46) 


1 


If  the  incident  electric  field  vector  E is  chosen  so  as  to  have  only  a 
positive  y-eomponent,  then  the  incident  electric  and  magnetic  field  vec- 
tors (E  , H*,)  can  be  expressed  as 


—I  -ikz  , „ 

E = E„ e = l E 


X . A ,, 

.1-  0 


(1) 


“0  *y  “0  ~ "xO  1 0'^* 


(47) 


i 

I 
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and 


4 - Hce-ikZ  . i*  H„  £ •Vo!1)'y0(V?)' 

X-  V 

Taking  the  curl  of  Equations  (47)  and  of  (48)  yields 


(48) 


VXEy  - E0  ^ A/oVx(L(1Vio) 


and 


VXH  = H 
-x  0 


Jx  aAVxU5 


(49) 


(5  0) 


-iwt 


Assuming  a time  dependence  of  the  form  e . Maxwell's  equations  for 
free  space  become 


e!  = ^7  x g 
k x 


and 


I i ^ T 

H =tVxE 
~ x k ~y 

so  that  E*  and  are  expressible  as 


(51) 


/5?.\ 
\ * 


E1  = — E 5 A „ \7  * 1 V i X i 
-y  fc  ^0  f-*  i 0 V y.O  -y 


(53) 


ana 


(54) 


In  terms  of  a solution  M ^ for  the  vector  wave  equation.  Equations  (53) 
and  (54)  become 


By  Z.Aj2% 


.JC-  o 


teo 


(55) 


r t 

1 U 


I 

u. 
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H1  = 

“"X 


Using  the  relations  between  the  unit  vectors  of  the  oblate  spheroidal 
coordinate  system  and  the  Cartesian  system  (Eqs.  (3a),  (3b),  and  (3c)), 

the  components  of  the  adenoidal  vector  become 

C,V*>  ■*’%,**. 

- "ifj  F'1^  + ? </  * f2)'!/2  sin»  Ui0  0,)  (,) 

+ t?F  - n"  </ ’ s2)  1/22  sin*  (1>V^0  («)  Ilf) 

X 

+ L F"1  ['<*  - /><>  ' + ?2>'‘  cos*  <UVo  (i>  <^D/o  <1> 

(1  + S2)<  if  + iV*  COS4  U ,y'n  fc,i  i(1,V /n  (S)  | (57) 


(V/n 

-1  \f.  ~~?L  , 2 2. -1/2 

= I„  F 1/1+9  (>?  + $ J 


cos^  u/rt  (n)  -rr*1^  (?) 

v \ u/  * U 


</  + *r1/‘cos*'1,VJ0(f)  ^U/OT 

+ '{  ' * * # sU"1’  '(■'0l?,dh70(,ll) 


2 2 -1/2 


^(1  +/)(/  + I2)"1  sin4U/0U>^(1,V/0<S> 
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Renee,  the  incident  electric  and  magnetic  field  (E*  , H* ) can  now  be  ex- 
pressed explicitly  as 


r 


/,  2 2 ,2-1/2  , (1)„  , , d „ 

A/o[_r  A ^ + * * sin,,i  V/n  ^ dhUYo  V 


+ is  / . A 

’'JPo  — - 


d 


V «£>*[*»  ‘ 'fM’f + «2,‘l  co»*  (1>vio «» jruio  'V 

v = o l 


2 2-1 


d (1), 


+1(l+f  >(,,  + f ) cos*U/o(^)^-'  'V/0  (?>jj 


(59) 


and 


r'  co 

<=  4Eo  P_1  Vo  ft*?  V * f2,"1/Z  COE  * V o <7»  4!1>V  o «5> 

OU  

5T  . l/  2 2 , 2 -1/2  (1>„  . . d 

-!|Z  Ai,M  (\  + f > cos  * Vo  <5  > dTi  Vo  v 

■ £ = 0 N *N 

OQ 

+ 4 27  Aio  [*>  - 1 V + f2*'1  ata*  (1>Vo  «f»  ^ Vo  <«> 

•N 

+ ^(1  + f){^  + ^)'A  sin^^0(^)^(1Hvo{?)  j-,  (60} 
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THE  SCATTERED  ELECTRIC  FIELD 


The  expression  for  the  scattered  electric  field  must  satisfy 
Maxwell's  equations,  the  boundary  conditions  at  the  surface  of  the  ohiate 
spheroid,  and  the  boundary  conditions  at  infinity.  In  order  to  satisfy 
the  Maxwell  equations  it  is  necessary  that  the  electric  field  vector  sat- 
isfy the  vector  wave  equation  and  the  diver  gence  condition.  At  the  sur- 
face of  the  scatterer  the  tangential  component  of  the  total  elect  r ic;  field 
must  vanish,  and  at  large  distances  from  the  scatterer  the  scattered 
electric  field  must  have  the  form  of  a spherically  divergent  wave. 


Inspection  of  the  asymptotic  form  of  the  radial  function  of  the  third 

(3)  (a)  {2> 

kind,  v ‘ , shows  that  this  combination  of  V,  and  ‘VL  meets 

./m  /m 

the  diverging  wave  condition.  To  satisfy  the  boundary  conditions  at  the 

surface  of  the  scatterer  as  simply  as  possible  the  values  of  m are 

chosen  so  that  the  scattered  wave  depends  on  <t>  in  the  same  way  as  the 

incident  wave. 

The  scattered  electric  field  is  then  expressed  explicitly  as 


_s  Ep  / (3)  x , . (3)  z \ 

- =i (v  ~/0  + b -/I  ) 


(3)  x 13)  z 

where  the  arK*  ’ — i are  inverse^>'  proportional  to  length, 

k = 2tt/\,  and  the  and  p ^ are  arbitrary  constants  to  be  determined  by 

the  boundary  conditions  at  the  surface  of  the  scatterer.  If  the  super- 
script (i)  is  replaced  by  the  superscript  (3)  in  Equation  (57)  the  expres- 
(3)  x 

sion  ror  Mpp  1S 
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«3L,x 


*/0 


■l„  f'!  \/l  + {2(,2  + t2)-i/Z6iniU/0(n)5i(3,V^0«) 


i F*1  l/l  -Y,  2 (nc  + cZf!//2  d *3  *V . (f)  — U„  fnl 

“?  r i ' ’ ’ ' JO  ' d'l  XQ  ' 


-if  2 2 2-1  fat  u 

i.  F !Ml-*>(n  + t)  cos  ^ ' ‘VjQ  (U~-  U/Q  (n) 


-a  ~ j 


2 2 

. - / 1 i • - \ i,. 

I •(  I » ! H / • '( 


2 -1 


I £ ) 


IJ  / r. ) — - ^ 3 It  \ 

JO  ' - ' cU  Vo"' 


(62) 


Also. 


O'  z _ . „1 

— /l  -n  F 


2.-1/2  2 . 2-1/2 


(3), 


Ml-O  (n  + t sin  * Vjl  (n)  '"V,^  (O 


+ i 


- ' 1 2-1/2  2 2 - 1 /■}  i \ \ 

L F -UiM  ) 7 ("  + O " sir,  * U,.,  (M'  V„,  U) 

s l t l t I 


+ i*  F 


Mi  - *Z)l/Z  (1  + i.Z)l/Z^Z  + O'1  cos  /3)V^  <« )~~  UJl  (u) 


? 1 / ? 2 1-6  2 2-1  ri  / -i  \ 

5 (1  - n ) '"(l  + * ) '“(  n + t ) cos^U.jtTO—^V^jU) 


(63) 


I S 

Since  the  total  electric  vector  at  any  point  is  given  by  E_  + E , the 
boundary  conditions  at  the  surface  can  be  expressed  as 


i S 
, E.  + E. 

I i . i 


= 0 


(64) 


L -H  = S, 


f 


i 
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r„i  „s  i 

E.  + E.  I = 0 

4 4J 

L 4=4 


Substituting  the  expressions  for  the  incident  and  scattered  fields 
(Eqs  59  and  61 ) on  the  surface  and  using  the  above  expressions  for 

O'  x (3)  z 

']V|^q  and  Equations  (64)  and  (65)  become 

‘ £ *Jo  " + so,1/2  ^o‘*>dT“,Y>o«»o> 


• £ {-,0  + 4,l/2  'Vc  «’> 5T  ‘\o  “o> 


/)-  n v 


+ id  - 1^1  (n>  <3>v^i  !v} 


Ao[ 

d-  n ^ u 


0 


+ " <!  + 5o>  u/o,n,dT(1)vo  ! 


lt  ) 

’o' 


j£{ 


V so  !1  - ”Z)  <3>vZo  ‘VdV 


+ t <1  + 5o>  %fO;rr'5!uo(io)  I 


+ + <3) 
-C  )\  u 


vn'VdT  (0 


i 
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fo (1  (‘  * fo)I/2  (3,V/><V 


(6?) 


The  above  system  of  equations  can  be  solved  for  a „ and  (V  by  use  of  the 

i A- 

orthop-nnal  n^onerties  of  Ihe  an.pular  functions  Ui  (/? ).  If  Equations 
° “ Jc  m ' I ' 

(66)  and  (67)  are  multiplied  by  and  integrated  over  the  interval 

- 1 < < i , then 


1 JL  A/0  (i  * { 
£-  0 


2,ia  r _d 

c [ d s 


in,.  ,,  . 

VW 


o»i/2i 

-v  V 


a (3),r  . . I , (3),.  , . , , I 

? v/o<V!  !i  + ^ ve i^o)1zf 


J 

(68) 


and 


> a - /n 

*6—'  70  I * ’ 

.€=  0 k 


r » 1 1 ^ 

f Si  _y_  **  >v  ( fc  ) i i + $ VA  V.  fjf  ) t y 

’ 0 ' I d £ U0  V 1 0 l 0 ’V  '4  f 

U-  -J  s 


= A |q  - hi  4£2)f  {3)V,  («  ))  I + Sr(3)V?0  (f  ) I 

i u \ u ? c 0 0 / 3 0 x 0 0 -j 


£=  0 


<l  +.’0,1/2  (df  '”V2-l(SoO  '5  - 11  ♦*;>‘/i,3V,  '?on6 


}’ 


c\  \ 

7 / 


where 


*i  W-lVjc  ('T)UL0 

f 1 ^ - 1 / 2 

Tz  (/>L)  J jV’1  "A  ' U/0  (r  UL0  {>0dX 
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x3  iv*‘  L)  = j_  x7  U/0  ' UL0  7 ' d:( 

ViL|!0'  -I^Ulo  «>4Uio 

«5  '‘■W-i,*1 


h U 


L L)  - f 


11  -V>,/2^0  !'( ! dTT  u(  1 ^>d'l  . 


The  six  integrais  above  can  be  evaluated  m terms  of  the  Beta 
functions  B (t,k).  These  functions  are  defined  by: 


1 


,2t..  2 ,k 


2k+i  k! 


(1  - ^ ) U,|  - B{t  + 1/2, k+1)  - -2t  + , J <27  4-  3)  . . { 2k  + 2t  + 1 ) 

where  k and  t are  non- negative  integers. 

For  even,  L even. 


( — r.  \ 

I u ; 


*'  * 1 


\ CO  oo 

* V “• 


k - 0 s = 0 


oo  tXJ 


= 0 s = 0 


,1 

r*  co 

f — » 

l 

> 

" J-l 

_ k = 0 

C2k 

c“ 

o - 

,2)k 

Zk  ■■ 

c,S 

,L0 

f!„ 

Z,k 
- r;  ) 

l 

2K  - 2s  J-l 

<1  <i  -M  J c“  li -/)k 

1 J L k - 0 “K 


dr? 

t 


dri 


*E  and  O refer  to  even  and  odd  values  respectively  cf  -/  or  L.. 
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=c 

■ Z 


y ' cf ° -lo 


k = 0 s = 0 


OO  IV 


2s  2k -2s 


y y c*>  CL0 

k = 0 s = 0 


B(l/2,k+l) 


k+1  , . 

2 k: 


4k  - 2s  3-  5 • • (2k  -i-  I ) 


(71  a.) 


For  0.  odd,  L odd, 


I,  (i,  L)  = I,  (O,  O) 

1 X 


_ i1r\  V fy 

V 

, . _ 2s  2k  - 2s 

k = 0 s = 0 


= ^ ^ 


V' 


XV 

'Zc 


x'O  rI,0 


13  (3/2,  k+1) 


2k+1  k! 


k = 0 s = 0 


2s  ~ 2k  - 2s  3-  5 • • • (2k  + 3) 


(?lb) 


Since  the  angular  eigenfunctions,  U,^  (/{ ),  are  orthogonal  in  the 
interval,  (-1,1),  ^ *** 


J1  ^L)=  */L 


(71c) 


where  6 u is  the  Kronecker  delta  function  and  N is  given  in  Equations 

,-L  xj  x-  w 

(28)  ar.d  (29). 

Using  a technique  similar  to  that  in  the  evaluation  of  I the  remain- 
ing integrals  become: 


oo  k 

V1 


I (E.B>.I  (0.0).  2*  Z Czk-2S  B (3/2, k-rl)  ,72a) 

k = 0 s ~ 0 


I.  (E,  O)  = I„  (O,  E)  = 0; 

£.  c 


(72b) 
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► 

? 

I I3(E,E)  = 

1 

I3  (E,  O)  = 

T IP  FI: 

*4  " 

I.  (E,  O)  = 
4 

I4(O.E)  = 
I5  (E,  E)  = 
I.  (E,  O)  = 

15  (O.E)  = 

16  (E.E)  = 

1/  (E,  O)  = 

o 

! i 

» 

i 
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I3  (0,0)  = 0 


1,(0, E)>  27  27^7  ,3  B (1/2,  k+l); 

k = 0 s = 0 


r m ni  = n 

-4  ■ - • - • 


r y . ,j° 


s C ~ C R (7.  /•>  L-  J.n 

w2s  2k  - 2s  ^ 


2 X X»*l  -i',CtU  C™ 


k ~ 0 s = 0 


->\r  - y 


4S 


P (3/2,  k+l)  ; 


i5  (O.  O)  = 0 


*0  rv 

Z Z< 


..LO 


k = 0 s - 0 
x k 


2 s «*k  - 2 s 


B (3/2, k+l) 


\ V £0  LO  . 

/ C r bo/?.. 

' 2s  2k  - 2,s  ' ' 


k = 0 s = 0 


k+l}; 


I&  (O,  O;  = 0 


OG 

Z 


2L  C2s  C2k-2s  i.S(3/2,k+2)-(2s+l)B(5/2, 
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(73a) 

(73b) 

(74a) 

(74b) 
(74c) 
<7  5a) 
(75b) 
(75c) 
(76a) 

k+l)] 

(76b) 


I 


WILLOW  RUN  RESEARCH  CENTER  - UNIVERSITY  OF  MICHIGAN 


UMM  116 


G (E,  E)  - I (O,  O)  = 0 
.>  3 


(73a) 


op  k 

V"  i'Q  LO 

I3(E.0).I  (O.E).  2,  2,  <=2.  C2k  -2s 

k = 0 S = 0 


B'l/2,k+l);  (73b) 


t /tt>  vr  \ - T ir\  m - n 
• t uji  '-‘i  x'-'i  ^ r u 

“±  -i 

co  k 

I (E,  O)  = -2  AT  s cf°  CLU  . B (?)  /z.  k+i) 
k = 0 s = U 


(74a) 


(74b) 


I (O.E) -2  A A (k  + 1 -S)  c{°  <A° 

4 <iK  ■ 

k = 0 s = 0 


B (3/2tk+l)  ; (74c) 


I.  (E,  E)  = I=  (G,  O)  = 


(75a) 


I_  (E.  O)  = 

3 


X C21c  - 2s  B<V2.k+D 


k = o s = o 


(75b) 


I5  (O.E)  = 


co  rw 

\ T V ^ /)r\  t n 

A /.  ct”  B(i/2 . k+1) : 

^ 2s  2k  - 2s  ' ‘ 


(75c) 


E (E.  E ) - I6  (O.O)  - 0 


(76a) 


I6  (E.O)  = 


y-  s ^ 6 1 L,o  r i 

''  C2s  ^2k-2s  L5^3/2, ,K+2)“(2s+l)B(5/2,k-ll)  | 

k - 0 s ::  0 3 

(76b) 


i 
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7 (O  - 


CO 

Z 


k 

V 


(2s+i  ) B (3/2 . k+i) 


k = 0 s = 0 


To  simplify  the  writing  cf  Equations  (68)  and  (69),  let 


(76c) 


t,  ■>  i /?  r a \ 

B.  = A,r  (1  +$:)-' - i -t:  ' 'v  «*n)  i.  (i,  l) 
-I  -eu  u j d%  0 0 1 


^-<1+4‘1/2  [sr  (3>v/o^o>]  !i 


? / ->  \ 
u . - 


'V,  <»„>  >,  V-  L> 


* A^C  " ^0!‘/2  [di  (‘Vo  «0>]  *J  + ll>Vi» 

1 + # o>  i if  ‘'So  <?0>J  6 *'■  L!  + '0  '^0  <V  '4  L> 

L , ,,  ,2,1/2  r d (3 J ..  .]  _ . , . . 

\£  ~ ~ *0  ''l  r ' G*  | _dj  Vo  '^O^J  *5  U'L,, 


Mi+j;)1/z  (3,V,0(#0)I6(AL). 


By  use  of  this  notation.  Equations  (68)  and  (69)  take  the  form 

!»,  Djr  + Oj)  • 2, 


0 


at, 

ZL  <a/  + P,  T/>  - (L  = 0,  1:  2 . . ) 

/y_r»  ■*  ' /i  n *-■' 


V. 

*? 


¥ 
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To  obtain  an  approximation  to  the  solution  of  these  tv/o  independent 
sets  of  equations,  write  the  first  four  equations  of  each  set  as: 


an  D +0  + 0 + 0 

U 0 


0 0 0 
. + ?0  G0  + 0 + P2  G2  + 0 + • • • = B0 


0 + Dj  + 0 + 0 + ,.+  0 + pj  Gj  + 0 -r  - . 


2 2 „ 2 
0 + 0 + a.  D_  + 0 + . . . + p _ G„  + 0 + p.  G_  + 0 . 
L Z u u /,  c 


0 + 0 + 0 + a,  D„  + 0 + . . 
3 3 


_3  _ „3 

. + p + u + p3  u + 


B 


= B_ 


. = b: 


and 


.0  0 0 0 
0 + a1  bj  + 0 + + 0 + . . .+  0 +Pj  F?  + 0 + 


o + 0 + a,  S,  + 0 + . . 
U 0 2 2 


2 2 

0 + a,  S.  + 0 + a S + 0 + . . 
11  3 3 


an  Sn  + 0 + Sl  + 0 + ' 


0 0 

= 0 + R +0  + R + .. 
l 3 


+ P0  T0  + 0 + T2  + 0 


= R1  + 0 + R1  + 0 + 
0 2 


0 + p T“  + 0 + p T + 0 + 


2 2 

= 0 + Rj  + 0 + + . . . 


3 3 

+ 0 + 0.  T.  t 0 + . 
K0  0 2 


= R'  + 0 + R + 0 + 


(SO) 


£ 
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Here  a,  and  a (n  = 0,  1,  . . . ) occur  in  the  first,  third,  . . . 

2 n 2n 

equations  of  system  (80)  and  in  the  second,  fourth,  . . equations  of 
system  (81).  In  particular,  as  a first  approximation,  to  ctn  and  (3^, 

consider  only  the  first  equations  containing  and  in  each  set. 

Then 


I 0 

! Bo 

i o1 

i “0 

>°o 


0 I 0 

i D, 


0 ! 

,i  j 

_0_  i 

,0  i 
7n  i 


and 


<V 


i he  a_  , and 
2n+l 


P,, 


equations  of  system  (80)  and  in  the  first,  third,  . . . equations  of  system 
(81).  Hence,  by  a similar  process  for  n = 0, 


si 

0 

To 

(82) 

! n 
!D° 

< 

1 .1 

^ a 

lao 

0 

Dr 

V 

g° 

* 

< 

(83) 

(n  = 0, 

* > 1 ' 

. ) occur  in  the  second,  fourth,  . . . 

1 

I 

B, 

G, 

1 

1 

n 

0 

r: 

T, 

i 

i 

! i 

1 

D, 

G, 

! 1 

1 

i o 

0 

S , 

T 

1 

(84) 
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For  a general  a or  p (n  = 1,  i,  . . . ) ihe  n + 1 equation?  from 

bii  £ n 

each  system  in  which  a,  and  0.  appear  must  be  used.  Similarly,  for 

In  cn 

a general  a_  , or  p,  (n  = i , 2,  . . . ) n + I equations  from  esch  system 
2nl-l  tn+i 

in  which  a'  , and  p.  , appear  must  be  used. 

Zn-fl  Zn+1 

Ij  Li 

Thus  in  the  above  expressions  for  a and  p the  quantities  B D... , 
LL.L.L  nn  x-  x 

G j,  R^;,  S^,  and  Tj  each  depend  upon  the  oblate  spheroidal  coefficients 

C^.  These  spheroidal  coefficients  are  completely  determined  by  the 

2 K * . . 

7/Tn 

defining  equations  for  C~  and  and  the  recursion  formulas  (I6)and 

(17). 


„s 

Ed 


jo  y ('  <3)„x 

k ^,o 

x-  0 


tA  (3'^Z  ^ 

' v<?  - >1  y 


fhi  \ 


where  '3,MX,  if.  given  by  Equation  (62)  and  '3,MZ  by  Equation  (63), 

k ^ J\ 


Substituting  the  asymptotic  forms  of  ,3  *V 


/m’  I 


Eqs.  (36)  and  (37) 


1 


into  Equations  (62)  and  (63),  and  letting  q = 1,  leads  to 

ikr 


and 


0)M* 

-JO 


{3)MZ.  = 2* 


1 e 

Ym^V°s  ^Jqio  T 


1 k 1 [ -i,?  sin  4 + i.^  cos  *] 


(86) 


ikr 


(87) 


For  q - i,  equations  (3a)  and  (3c)  reduce  to 


i_  = -i  cos  <i>  - i sin  * 


*'( 


~x 


(88) 


ana 


i = -i  sin  <f>  +•  i cos  <£  . 
"x  -y 


(89) 
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V 

ASYMPTOTIC  FORM  OF  THE  SCATTERED  ELECTRIC  F IELD 

In  the  previous  chapter  an  expression  was  obtained  for  the  scattered 
electromagnetic  field  everywhere  in  space,  The  back -scattering  cross- 
section,  however,  depends  explicitly  only  on  the  nature  of  the  electro- 
magnetic field  in  the  vicinity  of  the  source,  that  is,  in  the  asymptotic 
behavior  of  the  scattered  field  as  e -V  °°  and  >f  = 1 . 

<: 

The  expansion  of  E was  given  in  the  form 
E 
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So  that,  iiiquaiioris  (86)  and  '87  > become 

ikr 

<3,M*„  » i a i 

~Jl  0 'y  ‘iO  r 


(90) 


l/*i 


-1  -1  e 


ikr 


M„5*2ivq/i«  k 


(91) 


Substituting  Equations  (90)  and  (91)  into  Equation  (61)  and  neglecting 
. 2 

term  of  order  i/r  the  asymptotic  form  of  the  scattered  electric  field 
is  giv^n  by 


Ert  ikr 
£s  _ . 0^  e 

— ~y  k r 


V" 

JhO 


*1  qi0 


(92) 


-l.  I 
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VI 


TIIE  BACK -SCATTERING  CROSS-SECTION 
OF  THE  OBLATE  SPHEKOID 

The  asymptotic  expression  for  the  scattered  electric  field  (Eq.  92) 
makes  it  possible  to  formulate  the  far  zone  back-scattering  cross- 
section  explicitly. 

The  scattering  cross-section,  cr,  is  defined  as  the  cross-section  of 
an  isotropic  scatterer  which  would  scatter  in  the  direction  of  interest 
the  same  amount  of  energy  as  the  oblate  spheroid  scatters  in  that  direc- 

T 

tion.  If  P1  denotes  *he  total  power  intercepted  by  an  isotropic  scatterer 
of  cross-section  cr,  then 


where  p is  the  Poynting  vector  of  the  incident,  wave. 


(93) 


Thus  the  magnitude  of  the  Poynting  vector  for  the  wave  scattered  by 
the  isotropic  scatterer  is 


I 


4 -nr 


!~i  E 

1/0  0 

4rr  V u 2 " 

“ U3  r 


(94) 


where  r is  the  distance  from  the  center  of  the  isotropic  scatterer  to  the 
point  of  observation. 

But  j J is  also  given  by 


s 

s 

s 

p 

E 

X H 

V 


E 


V 


rrr  E„ 
/ u 0 


^0  k2r2!  £-  0 


> 


Q/9  9 /)(\ 
jc.  u 


(95) 
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So  that,  combir.ing  Equations  (94)  and  (95),  <r  is  expressed  in  terms  of 

the  constants  u-  of  the  scattered  electric  field  as 
x 


-*‘(T 

T 


or 


oo 

\ ’ 1 

_ ,2  | 

^ ! 

6 

JmC’*  2 i)  ^ /If-.  \ 

j,n  x \ 

7 r 

^ Q,;  Vn  1 
£ --0  i 

t 

(v  6 ) 

<r 

4vFZ 

~z~ 

f- 

— ^ 

Z 

(96a) 

where  q.  is  given  in  Equations  (34)  and  (35),  F is  one  half  the  focal 
^ ?-rrF 

length  and  £ - ~ — , X being  the  wavelength  of  the  incident  radia-oon. 

\ 

Leitner  and  Spence  (Ref.  3)  have  tabulated  q for  f = i,  Z,  3,  4, 

jl  0 

and  5 and  for  ^=0,1 , ...  5.  This  tabulation  of  a , must  be  extended 

-c.  ” 

to  larger  values  oiJ  and  t if  numerical  values  of  <r  are  to  be  obtained. 
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APPENDIX  I 

SCATTERING  FROM  A CIRCULAR  DISK 

j 

In  the  limit,  5A->0,  the  oblate  spheroid  approaches  a circular  disk 

\j 

of  radius  F.  Bouwkamp  (Ref.  9)  has  shown  that  the  scattered  field  from 
such  a disk  will  be  infinite  at  its  edge.  Flammer  (Ref.  19)  has  solved 
the  disk  problem  in  terms  of  oblate  spheroidal  wave  functions  which 
possess  the  correct  singularities.  To  demonstrate  that  the  scattered 
wave  has  the  correct  singularity  at  U = 0,  the  expansion  of  the  radial 
functions  in  terms  of  Hankel  functions  used  in  this  paper  is  inappropriate 
and  Equations  (34)  and  (46)  of  Reference  3,  giving  the  functions  as  power 
series  in£,  should  be  used. 


Then  it  can  be  shown  that  Equation  (b)  of  Part  II  in  Reference  10  and 
Equation  (6 1 ) of  this  report  differ  essentially  only  in  the  presence  of 
even  values  of-^in  the  summation  of  Equation  (61).  However,  even 

# 3 ) x 

values  of  / imply  even  values  of/-m  in  ' 'M^  and  odd  values  of-/-m  in  j 

(3)  z 

M ,, . It  can  be  shown  that  the  radial  functions  of  the  first  kind  vanish 
ati  = 0,  for  J-xrv  odd,  while  their  derivatives  vanish  at  § = 0 for-^-m 
even.  This  implies  that  B « = = 0 for  ■£  even,  in  Equation  (?7).  Thus 

Equations  (80)  for  the  and  with  J even  are  homogeneous  equations 
in  the  disk  limit,  with  the  solutions 


i 

) 


= 0 . 


Hence  (61 ) does  in  fact  reduce  to  an  odd  sum  ovevJ?.  The  singularities 
are  discussed  in  detail  in  Reference  10. 


The  edge  condition  of  vanishing  tangential  electric  field,  which  was 

used  by  Flammer  to  obtain  a unique  solution,  follows  automatically 

from  the  solution  presented  here  which  required  - 0 for  all  values 

of  . 
u 
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.APPENDIX  II 

THE  SCAT TERTNC-  CROSS- SECTION 
OF  THE  PROLATE  SPHEROID 

The  oblate  spheroid  can  be  transformed  into  the  prolate  spheroid 
by  the  simple  transformation: 

oblate  - prolate, 

§ oblate  = i % prolate, 

<!•  ohlate  = ir-#  prolate,  and 

e oblate  = i«  prolate.  (A) 


Therefore,  if  the  oblate  spheroidal  functions  and  oblate  spheroidal 
coefficients  are  transformed  into  prolate  spheroidal  functions  and  pro- 
late spheroidal  coefficients  respectively,  the  solution  of  the  boundary 
value  problem  for  the  oblate  spheroid  should  be  transformed  into  the 
solution  of  the  same  boundary  value  problem  for  the  prolate  spheroid. 


This  transformation  of  solutions  does,  in  fact,  occur.  When  the 
coordinate  transformation  is  applied,  the  oblate  spheroid  becomes  a 
prolate  spheroid  and  the  separated  angular  and  radial  differential  equa- 
tions are  arrived  at  from  Equation  (12).  These  equations  are  satisfied 
by  the  prolate  spheroidal  functions.  Furthermore,  the  recursion  formu- 
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las  for  the  coefficients  C, , take  the  form, 
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Because  these  recursion  formulas  differ  for  the  oblate  and  prolate 
spheroidal  systems,  the  spheroidal  coefficients  are  different,  numerically, 
in  the  two  systems.  This  numerical  difference  leads  to  different  values 
of  the  scattered  electromagnetic  field  and  of  the  scattering  cross-section 
even  though  the  explicit  form  of  these  two  quantities  is  the  same  in  both 
coordinate  systems. 
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